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ABSTRACT 



The chemical properties of high-z galaxies provide important information for constraining galaxy evolutionary scenarios. However, 
widely used metallicity diagnostics based on rest-frame optical emission lines are unusable for heavily dust-enshrouded galaxies (such 
as submillimeter galaxies; SMGs), especially at z > 3. Here we focus on the flux ratio of the far-infrared fine-structure emission lines 
[N n] 205 //m and [C n] 158 //m to assess the metallicity of high-z SMGs. Through ALMA cycle observations, we have detected the 
[N n] 205 fim emission in a strongly [C n]-emitting SMG, LESS J033229.4-275619 at z = 4.76. The velocity-integrated [N n]/[C n] 
flux ratio is 0.043 ± 0.008. This is the first measurement of the [N n]/[C n] flux ratio in high-z galaxies, and the inferred flux ratio 
is similar to the ratio observed in the nearby universe (~ 0.02 - 0.07). The velocity-integrated flux ratio and photoionization models 
suggest that the metallicity in this SMG is consistent with solar, implying that the chemical evolution has progressed very rapidly 
in this system at z = 4.76. We also obtain a tight upper limit on the CO(12-ll) transition, which translates into CO(12-l l)/CO(2-l) 
< 3.8 (3tr). This suggests that the molecular gas clouds in LESS J033229.4— 275619 are not significantly affected by the radiation 
field emitted by the AGN in this system. 

Key words, galaxies: high-redshift - galaxies: individual (LESS J033229.4-275619) - submillimeter: galaxies - submillimeter: ISM 



1. Introduction 

Investigating the metal content in galaxies is a powerful diagnos- 
tic for testing galaxy evolutionary models, because the metal- 
licity is determined by the past star-formation history, gas in- 
flow, and outflow phenomena. The gas metallicit y in galaxies has 
been investigated in ga laxies up to z~3 dMaiolino et alj l2008t 
I iMannucci et alJ l2010[). However, rest-f rame optical metallic- 



[C n] 158 //m emission (v rest = 1900.539 GHz) is one of the 
strongest lines among the whole electromagnetic spectrum and 
is con sequently observed for some high-z objects even up to z ~ 
6-7 dMaiolino et all 120051: 1 Walter et all l2009atlVenemans et al.1 
|2012|) . This emission arises mainly in H n regions and pho- 
todissociation regions (PDRs), and the relative contribution from 
those two regi ons depends on the physical conditions of gas 
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ity diagnostics (e.g.jNa gao et al.1 [2006al) are not applicable for 
galaxies at z > 4, where the optical emission lines required to 
measure the metallicity are redshifted out of th e near-infrared at- 
mospheric window (see Mai olino et al.1 120081) . Moreover, high-z 
young galaxies with vigorous star formation are often obscured 
by dust and recognized as, e.g., ultraluminous infrared galaxies 
(ULIRGs) and submillimeter galaxies (SMGs). In these galax- 
ies, the rest-frame optical lines are significantly affected by the 
dust extinction and accordingly the optical metallicity diagnos- 
tics may lead to large systematic errors. For instance, it has been 
claimed that ULIRGs show systematically lower metallicities 
than expected from the mass-metallicity relation dRupke et alJ 
2008: lCaputi et alJ [2008); but this may simply be a consequence 
of the fact that in these heavily obscured systems the optical 
metallicity tracers o nly probe the outer, less enriched regions. 
ISantini et alJ (12010) reported that the high dust mass measured 
in these galaxies disagrees with the metallicity inferred from the 
optical metallicity diagnostics, suggesting that the optical spec- 
tral indices are likely unreliable in these dusty systems. 

Metallicity diagnostics exploiting far-i nfrared fine-s t ructur e 
emission lines are powerful alternatives (INagao et al.l 1201 lh . 
Here we focus on the two strongest fine-structure lines at 
^ rest > 150 /an, i.e., [C n] 158 fan and [N 11] 205 /mi. The 



clouds (e.g.. lAbel etaj] 120051). T he TN 11] 205 fim emission 
(v res t = 1461.132 GHz ;lBrownetai1l 1994 arises in H 11 regions 



(see lNagao et all 120111) . In H 11 regions, the [N n]/[C 11] flux ratio 
depends mostly on the N/C elemental abundance ratio, while it 
is relatively insensitive to other gas physical conditions such as 
the gas density, since their critical densities are similar (n cl - = 44 
cm" 3 and 46 cnr 3 for [N 11] 205 fim and [C 11] 158/rni at 8000 
K, respectively). Note that the N/C ratio is proportional to Z gas 
at least for log Z pas /Z > -0.5 as observed in Galactic H n re- 
gions (e.g. . Ivan Zee et all 1 1998b because nitr ogen is a secondary 
element (see also Hamann & Ferland 1999). Since the volume 
ratio of H 11 regions and PDRs also depends on cloud proper- 
ties such as the gas density and ionization structure (hence the 
ionization parameter), the observed [N n]/[C 11] flux ratio ac- 
cordingly depends on those parameters. Recently, [N 11] 205 /mi 
detections have been ob t ained for three strongly lense d galaxies 
at z~4-5 dDecarli et al.1 l2012t ICombes et all 120121) : however 
a [C 11] 158//m detection has not been reported in any of these 
galaxies. 

Here we show that the [N n]/[C 11] flux ratio is a good 
metallicity indicator based on photoionization models, and apply 
this new method to a high-z luminous SMG, LESS J033229.4- 
275619 (hereafter LESS J0332). We selected this SMG be- 
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Fig. 1. Velocity-integrated [N n] map of LESS J0332 after the 
continuum subtraction, with the spatial sampling of 0.2 arc- 
sec/pixel. The velocity range from -558 km s to -154 km 
s _1 (with respect to the Lya emission-line peak) is integrated. 
Contours at 3cr, 4<x, 5<t, 6<x, and 7<x levels are also given in the 
map. The shape of the synthesized beam is shown at the lower 
left-hand corner. 



cause this object is st arburst-dominated and s hows intense 
[C n] 158/mi emission dDe Breuck et al.l 1201 ll) . making this 
SMG a good target for the [N n] 205 /mi observation. In this 
Letter we report a clear [N n] 205 pm detection based on our 
ALMA cycle observation, and discuss the chemical property 
of a high-z SMG. Throughout this Letter, we adopt a cosmology 
with H = 7010ns- 1 Mpc" 1 , Q m = 0.27, and Q A = 0.73. 



2. Observations and results 

We observed LESS J0332 at z = 4.76, with the band 6 receiver 
on the Atacama Large Millimeter/submillimeter Array (ALMA) 
in the dual-polarization setup, as a cycle program. The obser- 
vation was carried out in three separated runs; on 201 1 October 
1, 2012 January 12, and 2012 January 13. Each observing run 
consisted of 1.2 hours and the total observing time is 3.6 hours, 
including overheads. The receiver was tuned to 253.96001 GHz 
to cover the redshifted [N n] 205 pm emission with the upper 
sideband, whose redshift is based on our previous [C n] 158 pm 
detection dDe Breuck et al.l 1201 lb . The lower sideband was used 
to cover the redshifted CO( 12-1 1) emission (v rest = 1381.995 
GHz). The correlator was used in the frequency domain mode 
with a bandwidth of 1875 MHz (488.28 kHz x 3840 channels). 
The observation was performed with 18 antennas in the com- 
pact configuration, but the data obtained with one antenna was 
flagged out in the last run due to its high system temperature. 
Callisto was also observed as a flux calibrator. The bandpass and 
phase were calibrated with J0522-364 and J0403-360, respec- 
tively. The atmospheric condition was PWV = 3.0-5.0 mm in 
the first run and PWV = 1.0-2.0 mm for the last two runs. 

The data were processed with Common As tr onomy 
Software Applications (CASA: iMcMullin et alj [2007; 
iPe try & the CAS A Development Team! 120121) in a standard 
manner. A 70-channel (or equivalently, 40.54 km s -1 ) binning 
was applied to the data cube, and then the clean process was ap- 
plied with the natural weighting, which gives a final synthesized 
beam size of 1.67" X 1.48" (position angle = 26.8 degree). The 
[N n] 205 pm emission of LESS J0332 is clearly detected in the 
continuum-subtracted binned channel map, as shown in Fig. 1 . 
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Fig. 2. Top panel: ATCA CO(2 -l) spectrum of LE SS J0332 
adopting a 90 km s _I binning (Cop pin et al.l l2Q10t) . Middle 
panel: APEX [C n] 158j um spectrum with a 28 km s _1 binning 
dDe Breuck et all 1201 ll) . Bottom panel: ALMA [N n] 205 pm 
spectrum with a 41 km s _1 binning. All the spectrum is 
continuum-subtracted. The best-fit Gaussian profile is overlaid 
on the [C n] and [N h] spectra. The velocity is relative to the Lya 
redshift, and the zero velocity is indicated by the dotted line. 



Fig. 2 shows the [N n] 205 pm spectrum of LESS J0332 
adopting an aperture size of 4.8 arc sec 2 , with the previously re- 
ported spectra of CO(2-l ) (ICoppin et al.ll2010b and [C n] 158 pm 
De Breuck et all 1201 ll) . The [N n] 205 pm detection signifi- 
cance is ~8<t. The emission line is fitted with a single Gaussian 
profile. The best-fit result is parameterized by the peak fre- 
quency of 253.868+0.007 GHz, peak intensity of 2.58 + 0.21 
mJy, FWHM of 230 + 22 km s _1 , and the velocity-integrated 
flux of 630 + 78 mJy km s _1 . The redshift inferred from the ob- 
served [N n] 205 pm frequency is Z[mi] = 4.7555+0.0002, which 
is significantly blueshifted from the Lyor redshift. This blueshift 
is probably caused by the Lya absorption by the intergalactic 
neutral hydrogen. The redshift is consistent with the CO(2-l) 
redshift, though there is a velocity offset, significant at ~2<x, be- 
tween the [N n] 205 pm and [C n] 158 pm lines. This velocity 
offset may be a consequence of the modest signal-to-noise ratio 
in the two spectra, especially in the [C n] spectrum. Forthcoming 
[C n] ALMA observations of the same objects will either con- 
firm or reconcile the velocity discrepancy. However, if this ve- 
locity offset is confirmed, this results in different [N n]/[C n] flux 
ratios in the blue and red components (see §3.2). 

The redshifted CO(12-ll) emission is not detected in our 
ALMA data. The inferred 3<x upper limit on the velocity- 
integrated flux is 344 mJy km s _1 (adopting the same aperture 
size as adopted for the [N n] 205 pm line and assuming the same 
velocity profile). 

The continuum emission at A t, s ~ 1.2 mm (i.e., A rest ~ 
210 pm) is clearly detected with high accuracy. Its flux is 3.5 + 
0.1 mJy with the same aperture as adopted for the [N n] 205 pm 
measurement. These measurements were made by combining 
the upper and lower sideband data and excluding the channels 
affected by the [N n] 205 pm emission and by a relatively high 
noise level at the edges of the sidebands. 
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Table 1. Observed properties of LESS J033229.4-275619 



Parameter 


Value 


Reference 


ZLyo- 


4.762±0.002 


CooDin et al. (2009) 




4.755±0.001 


Coooin et al. (2010) 


ztcrrj 


4.7534+0.0009 


De Breuck et al. (2011) 


Z[N1I] 


4. /333+0.0U02 


this work 


A:o(2-i) 


0.09+0.02 Jy km s -1 


CoDDin etal. (2010) 


/[crrj 


14.7±2.2 Jy km $~ l 


De Breuck etal. (2011) 




0.630±0.078 Jy km s" 1 


this work 


/C0(12-ll) 


< 0.344 Jy km s" 1 (3<r) 


this paper 


AVcop-i) 


160±65 km s _I 


CoDDin et al. (2010) 


AV[ciij 


161±45kms-' 


De Breuck etal. (2011) 


AV [N „] 


230±22 km s" 1 


this work 



Both the [N n] 205 fim emission and the continuum emission 
are unresolved at our angular resolution. The two-dimensional 
Gaussian fit on the velocity-integrated images results in the ma- 
jor and minor axis FWHMs of 1 .86 + 0.44 arcsec and 1 .3 1 ± 0.47 
arcsec for the [N n] 205 fim emission, and 1 .72 + 0. 19 arcsec and 
1.55 + 0.21 arcsec for the continuum emission. Both results are 
consistent with the synthesized-beam image shape. 



3. Discussion 

3.1. Possible AGN contribution 



iGilli et all d201 lb reported the presence of a Compton-thick ac- 
tive galactic nucleus (AGN) in LESS J0332. If the AGN con- 
tributes significantly to the observed far-IR emission lines, its 
interpretation would become accordingly more complex. Within 
this context the CO spectral line energy distribution can help, 
since it is sensitive to the heating energy source. CO lines at high 
excitation levels are significantly stronger when the molecular 
gas clouds are affected by the X-r ay emission from AGNs tha n 
in cases without AGNs (see, e.g.. ISpaans & Meijerinkl [2008). 
The nearby ULIRG-Quasar Mrk 231 shows strong high-/ CO 
lines up to J = 13-12, which are properly accounted for 
by introducing X-r a y do minated regions (XDRs) into models 
( Ivan der Werf et all l2010h . At high-z the gravitationally mag- 
nified quasar at z = 3.91, APM 08279+5255 also show s strong 
high-/ CO lines dWeifi et all l2007tlBradford et all 1201 ll) . which 
are also well described by XDR models. The CO spectral line 
energy distribution of other star-formation-domina ted high-z 
galax ies, such as SMM J16359+6612 at z = 2.5 (IWeiB et al.1 
l200l and IRAS F10214+4724 at z = 2.3 dAo et al l ^OoIk Ts 
completely different from that of APM 08 279+5255, showin g 
weaker high-/ CO lines (see, e.g., Fig. 14 in lWeiB et alJ 120071) . 

By combining our upper limit on the C O( 12-11) flux and th e 
previous measurement of the CO(2-l) flux (ICoppin et al.l 120101) . 
we obtain a 3<x upper limit on the flux ratio of CO(12-l l)/CO(2- 
1) of 3.8. This upper limit is inconsistent with the CO spectral 
line energy distribution of the quasar APM 08279+5255, but is 
fully consistent with other star-f o rmati on-dominated high-z ob- 
jects (see Fig. 14 in IWeiB et all 120071) . This suggests that the 
molecular clouds in LESS J0332 are not described by XDR mod- 
els, i.e., the AGN contribution to the heating and excitation of 
the ISM in LESS J0332 is not sign ificant. This is consiste nt with 
our earlier study on LESS J0332 dDe Breuck et al.l 1201 ll) , where 
we estimated that the XDR contribution to the [C 11] 158/mi is 
-1.3%, based on the absorptio n-corrected X-ray luminosity of 




^2-iokev = 2.5 x 10 44 erg s 1 (IGilli et al.l 12011 and a scaling 



relation of L[ciij,agn = 2 x 10 3 L2-iokeV dStacev etalJ 12010). 



z/z 

' sun 

Fig. 3. Observed [N n]/[C n] flux ratios compared with model 
predictions. The green hatched range denotes the observed range 
for low-z galaxies. The horizontal dashed line shows our ALMA 
result on LESS J0332, where the emission-line fluxes are mea- 
sured by integrating the best-fit Gaussian function. Dotted ma- 
genta lines denote the flux ratio at the red part (from -350 km 
s _1 to -150 km s _1 ) and blue part (from -550 km s _1 to -350 km 
s~') of the lines. The red and blue lines show Cloudy model re- 
sults as a function of Z gas with log «hii = 1-5 and 3.0 respectively, 
while solid and dashed lines denote the models with log t/nn = 
-2.5 and -3.5, respectively. 



3.2. Gas metallicity 

Based on our [N n] 205 fim detection an d our previous 
[C n] 158 fim detection dDe Breuck et all l201ll) in LESS J0332, 
the velocity-integrated flux ratio of [N n]/[C 11] is inferred to 
be 0.043 + 0.008. Unfortunately, there are only few previous 
measurements on the [N n] 205 fim line in galaxies (mostly be- 
cause this line was not covered by the ISO/LWS wavelength 
range). In the nearby univers e, the flux ratio of [N nl/f C 11] is re- 
ported o nly for M82 (~0.050 ; |Petuchowski et al.l 1 9941). Mrk 231 
(-0.0 67; iFischer et al. l|2010|), NGC 1097 (-0 . 017: iBeirao et ail 
[2010h . and Arp 220 (-0.059: iRaiigwala et al.ll201 ll) . Therefore 
the [N n]/[C n] flux ratio of LESS J0332 is similar to the ob- 
served ratios reported for nearby galaxies, suggesting similar 
7 

■^gas- 

At high-z [N n] 205 fim h as been dete c ted in HLS 
J091828.6+514223 at z = 5.24 dCombes et al.l 120121) . APM 
08279+5255 at z = 3.91, and MM 18423+5938 at z = 3.93 
dDecarli et al.l 120121) . However, their [C 11] 158 /mi line has 
not been observe d and consequently their [N n]/[C 11] ratio 
is unknown (see IWalter et al.l 12009b and references therein). 
Therefore our [N 11] 205 fim detection allows us to infer the first 
measurement of the diagnostic [N n]/[C 11] flux ratio at high- 
z. Note that th ere are many [N nl 122 / mi detections in nearby 
galaxies (e .g JGracia-Carpio et al.l 1201 ll) and also in a few high- 
z galaxies (Ferkinhoff et al] 1201 ll) . Although there are attempts 
to infer th e fN n] 205 fim flux from the [N 11] 122 fim emis- 
sion (e.g., iDecarli et al.1 120121) . this method may introduce a 
large systematic error because the flux ratio of [N n] 122 fim and 
[N n] 205 fim varies by a fact or of -10, i . e., it i s strongly depen- 
dent on the gas density (see lOberst et al.l 120061) . 
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To explore the gas metallicity in LESS J0332 more quan- 
titatively, we carrie d out model calculations using Cloudy 
(Ferl and et al.l 1998) version 08.00. Since the [C n] line arises 
in both H n regions and PDRs, a consistent treatment to connect 
those two regions is required to investigate the [N n]/[C n] flux 
ratio. We assumed a pressure-equilibrium gas cloud that is char- 
acterized by certain gas densities and ionization parameters at 
the illuminated face («hii and Uhu) for each model run. Here we 
examine gas clouds with log«Hii =1.5 and 3.0, and log Uhu = 
-3.5 and -2.5, for Z gas /Z = 0.05 - 3.0. Note that the gas density 
in PDRs is higher than log «hii under the pressure-equilibrium 
assumption. The input co ntinua are the st arburst99 constant star- 
formation model spectra dLeitherer et al.l Il999l) with an age of 1 
Myr and a stellar metallicity equal to the gas metallicity. We did 
not take AGN effects into account in the models since the ex- 
pected AGN contribution is small (see §3.1). The relative chem- 
ical composition of gas clouds is scaled from the solar elemen- 
tal abundance ratios except for nitrogen, w hich scales as Z gas at 
Z gas > 0.23Z o ("see lKewlev & Dopitall2002h . Orion-type graphite 
and silicate grains are included. Calculations are stopped at the 
depth of Ay = 100 mag to cover the who le [C n] 158/im emit- 
ting regions, following lAbel et all d2005l) . Note that the result- 
ing [N n]/[C n] flux ratios are not sensitive to the a dopted stop- 
ping c riterion. The calculation details are given in Naga o et al.l 
(1201 lh . 

The [N n]/[C n] flux ratio predicted by the model runs are 
compared with observations in Fig. 3. The models show that the 
[N n]/[C n] flux ratio increases monotonically with Z gas . The 
moderate dependences on parameters other than Z gas , such as 
«hii and Uhu, prevent us from determining the accurate Z gas 
value. The observed [N n]/[C n] flux ratio (~ 0.043) infers 
log(Z gas /Z ) ~ 0.0 + 0.3; i.e., consistent to the solar metallic- 
ity. Although the constraint on the gas metallicity is currently 
fairly loose, by observing additional lines in the future (such as 
[N n] 122 /mi and [O i] 145 //m) it will be possible to constrain 
the gas properties (density and ionization parameter) and there- 
fore determine the gas metallicity much more accurately. 

However, we note that the velocity profiles of the 
[N n] 205 fim and [C n] 158 /mi emission are different (Fig. 2). 
In the redder part of the line the [N n]/[C n] flux ratio is signifi- 
cantly lower than that in the bluer part. If the different [N n]/[C n] 
ratios for different velocities are caused by metallicity variations, 
this implies that the system is not chemically homogeneous, but 
probably resulting from a merging systems where gas with dif- 
ferent metallicities has not fully mixed yet. More specifically, the 
redder component with a high [N n]/[C n] flux ratio may be asso- 
ciated with a chemically enriched galaxy that is merging with a 
metal poor galaxy, associated with the bluer component charac- 
terized by the low [N n]/[C n] ratio. To investigate this scenario 
more quantitatively, we measured the [N n]/[C n] flux ratio on 
the redder part (-350 km s < V < -150 km s~') and on the 
bluer part (-550 km s < V < -350 km s~') of the lines, result- 
ing in a flux ratio of 0.086 and 0.020, respectively. These values 
roughly correspond to the highest and lowest flux ratios seen in 
nearby galaxies (Fig. 3), which is indicative of solar-supersolar 
metallicity in the former case and sub-solar metallicity in the 
later case. Note that the two distinct [N n] 205 /mi velocity com- 
ponents are seen in another high-z [N n] 205 /mi emitter, HLS 
J091828.6+514223 (Combes et al. 2012), which may be another 
example of a chemically inhomogeneous system at high-z. 

The low-metallicity component in LESS J0332 was already 
identified in our previous work on this SMG, where we ar- 
gued (based at that time solely on the [C n] 158 /mi, FIR and 
CO properties) that the gas metallicity in LESS J0332 is low 



(iDe Breuck et al.1 1201 lh . However, the most interesting result 
obtained here is that this system does also host a significant 
metal-rich component (with solar-supersolar metallicity), indi- 
cating that this system is chemically highly evolved already at 
z = 4.76 (the cosmic age of 1.27 Gyr). 

This result is consistent with past metallicity studies on high- 
z AGNs. The diagnostic flux ratios of some UV metallic per- 
mitted lines of type-1 quasars show no redshift evolution of 
the quasar broad-li ne region up to z ~ 6 (e.g., Nagao et al. 
2006cJ lJuarez etldl l2009t) . The metallicity in narrow-line re- 
gions of type-2 AGNs (which trace larger spatial scales and 
are therefore more closely related to the host galaxy properties 
than to the broad-line reg ion) also show no re dshift evolution 
(e.g.jNagao et alJ l2006bh even up to z ~ 5.2 (fMatsuok a et alJ 
120111) . Note that the high-z AGNs investigated i n those studies 
are ho sted mostly by massive g alaxies (see, e .g..|Sevmour et al.l 
120071) . similar to SMGs (e.g., iBorvs et al.l 120051) . These re- 
sults on the metallicity of AGNs and LESS J0332 suggest that 
the chemical evolution of massive systems has progressed very 
rapidly in the early epoch, z > 5, which is qualitatively consis- 
tent with the so-called downsizing evolution seen in g alaxies and 
AGNs (e.g.. lDeLuciaeTaII l2006lllkeda etaD l20ll . 
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